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Abstract. - We calculate the effective force mediated by thermally excited capillary waves 
between spherical or disklike colloids trapped at a fluid interface. This Casimir type interaction 
is shown to depend sensitively on the boundary conditions imposed at the three-phase contact 
line. For large distances between the colloids an unexpected cancellation of attractive and 
repulsive contributions is observed leading to a fluctuation force which decays algebraically 
very rapidly. For small separations the resulting force is rather strong and it may play an 
important role in two-dimensional colloid aggregation if direct van der Waals forces are weak. 



Introduction and synopsis. - The effective forces between rigid objects immersed in a 
fluctuating medium have attracted a steadily growing interest because their understanding 
allows one to design and tune them by choosing suitable media and boundary conditions 
and by varying the thermodynamic state of the medium. Possible applications range from 
micromechanical systems to colloidal suspensions and embedded biological macromolccules. 
Accordingly these fluctuations may be the zero-temperature, long-ranged quantum fluctu- 
ations of the electromagnetic fields giving rise to the original Casimir effect [1] between flat 
or corrugated immersed metallic bodies [2,3]. Other examples for fluctuation induced long- 
ranged effective forces between immersed objects involve media such as bulk fluids near their 
critical point [4], membranes [5] or interfaces [6]. 

In this work we investigate the latter manifestation of thermal Casimir forces for nanocol- 
loids floating at surface-tension dominated liquid-vapor or liquid-liquid interfaces where they 
are effectively trapped. In the presence of charges their mutual interactions often exhibit unex- 
pected strong and long-ranged attractions [7] inducing mesoscale pattern formation. However, 
these unusual patterns, which, e.g., are of interest for optical applications once they are fix- 
ated on a solid substrate, appear also for neutral nanocolloids [8]. In spite of some progress 
[9], the nature of these effective forces between the colloids is not yet fully understood. This 
pertains in particular to the role of fluctuations for these types of observations. 

Therefore the present effort aims at determining the fluctuation induced contribution to 
these forces generated by the inevitable thermally excited capillary waves of the fluid interface. 
As a first step the colloids are taken to be spherical or disklike and electroneutral. We shall 
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pay special attention to the fact that these colloids are of finite size and that, through their 
surface properties, they exert certain restrictions on the three-phase contact line formed at 
the intersection of the fluid interface with the colloid surface. 

For the idealized situation of a pinned contact line (corresponding to Dirichlet boundary 
conditions) this problem has already been studied, mostly in the context of membrane-type 
interfaces [10]; some general results can be inferred also from Ref. [11]. For a pinned contact 
line and a Gaussian weight for the interfacial height fluctuations these studies yield for the 
asymptotic decay of the fluctuation induced force F(d) as function of the center-to-center 
distance d between the colloids (i) F(d — > oo) oc dr 1 (attractive) if the centers of the colloids 
are fixed by some external means and (ii) F{d — > oo) oc d~ 5 (attractive) if the colloids are free 
to fluctuate vertically. Going beyond this analysis we find that F(d — ► oo) oc d~ 9 (attractive) 
if, corresponding to the generic situation, the fluctuations of the three-phase contact line are 
included, governed by a boundary Hamiltonian which is derived from the surface properties 
of the colloids and of the interface. This highlights the importance of boundary conditions in 
Casimir problems. In our case they render the long-ranged tail of the effective force virtually 
unimportant. However, if the colloids approach each other with a surface-to-surface distance 
h — » we find that the force increases oc h~ 3 / 2 . Thus for actual systems the fluctuation 
induced force is effectively short-ranged and therefore it may play an important role in the 
coagulation of interfacially trapped colloids. This behaviour cannot be captured if the colloid 
is approximated by an essentially pointlike object as it has been done in Ref. [6]. We shall 
also discuss experimental possibilities to separate the fluctuation induced forces from the 
ubiquitous dispersion forces acting in colloidal systems. 

Model. We consider an interface between two fluid phases I and II at which two 
nano- or microscopic colloids are trapped, either spherical with radius R or disklike with 
radius R and thickness H. Since the weight of such colloids is negligible, the equilibrium (or 
reference) configuration is the flat interface with the centers of the spherical colloids vertically 
positioned such that Young's law holds at the horizontal three-phase contact circle with radius 
ro = R sin 8. Young's angle is determined by cos 8 = (71 — 711V7 where 7 is the surface tension 
of the I/II interface, and Tipi] the surface tension between the colloid and phase I [II]. For 
the disks, the contact line is cither the upper (8 < 7r/2) or lower (8 > ir/2) circular edge, 
so that ro = R- Deviations from the planar reference interface z = arc assumed to be 
small which allows us to use the Monge representation (2, y, z = u(x, y)) = (x, z = u(x)) as a 
parametrization of the actual interface positions. The free energy cost for thermal fluctuations 
with small gradients around the flat interface is determined by the the change in interfacial 
energy of all interfaces (I/II, colloid/I [II]): 



n tot = ft™ + 71AA1 + 7uAA n 



? / 



■ 71AA1 + 711AA11 . (i) 



Here, 7i' C w is the standard capillary-wave Hamiltonian. In cq. the capillary length is given 
by A c = [7/dpn — pi\ g)] 1 ^ 2 , where pt is the mass density in phase k and g is the gravitational 
constant. The first term in W cw expresses the energy needed for creating the additional 
interface area associated with the height fluctuation, the second one the corresponding cost 
in gravitational energy Usually, one has A c 3> R] however, care is required in taking the limit 
A c — > 00 (see below), corresponding to the limit of a vanishing pinning force g — > 0. The 
integration domain S is given by the plane R 2 (z = 0) minus the enclosed areas <5i and S2 
of the projections of the contact lines on the first and second colloid, respectively, onto this 
plane. Thereby the Hamiltonian depends explicitly on the positions of the colloids. In the 
reference configuration, the projections Sf of the contact lines are the aformentioned circles 
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with boundary dS® and common radius r$. As discussed in the Introduction, in the case 
of spheres the contact line itself may fluctuate around its reference position with an energy 
cost which is determined by changes of the projected meniscus area S and of the interfaces 
colloid/1,11 (A.Ai,ii ^ in eq. Q) with respect to the reference configuration. We introduce 
the vertical position of the contact line at colloid i as a function of a polar angle ipi defined on 
dS® by fa = u(dSf) = P,o + Em>il^™ cosmcpi + Qi m s'm inept}. The expansion coefficients 
Pi m and Qim are referred to as multipole moments of the contact line. Furthermore, hi is 
introduced as the fluctuation induced change in vertical position of the colloid centers with 
respect to the reference configuration. Following Ref. [12], we restrict the integration in TC' CW 
to the reference domain R 2 \ (J Sf and expand the remaining energy differences up to second 
order in fa and h^, introducing a boundary Hamiltonian 7ib,i so that 7itot = H cw + J^,- 7ib,i' 

WcwM = n / d 2 x 

Wb,i[/i, K] = - f dip t (fi - hi) 2 = 7T7 
2 Jas° 

The effective force F(d) = — §j as function of the mean local distance between the colloid 
centers is determined by the free energy J-(d) = — ksT (In Z(d)). The partition function Z(d) 
is obtained by a functional integral over all possible interface configurations u and fa; the 
boundary configurations are included by S- function constraints, as introduced in Ref. [13]: 

Z = Z^ I p U cxp{-^i) fi / T>f< II SH Xi )-fa( Xi )]e X p\- Hh f^ hi] } . (4) 

Zq is a normalization factor such that Z(d — ► oo) = 1. We shall discuss two different realiza- 
tions of the boundary conditions for the contact line (see Fig. ^) . (A) The contact lines and the 
vertical colloid positions fluctuate freely; this corresponds to the physical situation of smooth, 
spherical colloids. In this case, the integration measure is given by Vfa = J dhi J7 m dPi m dQi m . 
(B) The contact lines are pinned in the reference configuration and do not fluctuate with re- 
spect to the colloid surfaces. This corresponds to disks or Janus spheres consisting of two 
different materials. Within the pinning case, we furthermore distinguish the following three 
physical situations: (Bl) The colloid positions are fixed, e.g., by using optical tweezers; thus 
there arc no integrations over the boundary terms. (B2) The vertical positions of the colloids 
fluctuate freely, so that Vfa = dhidPio5(Pio — hi). (B3) The vertical position and the orien- 
tation of the colloids (tilts) fluctuate freely. Up to second order in the tilts this corresponds 
to Ufa = dhidPiodPiidQn8(Pio — hi). The (^-function expresses the pinning condition. All 
these cases can be discussed conveniently after splitting the field u of the local interface po- 
sition into a mean-field and a fluctuation part, u = it m f + v. The mean-field part solves the 
Eulcr-Lagrange equation (—A + A f T 2 )u m f = with the boundary condition u m f \ggo = fa. 
Consequently the fluctuation part vanishes at the contact line: v \ 9S o = 0. Then the partition 
function Z = Za uc Z m f factorises - due to the Gaussian form of Ti. cw - into a product of a 
fluctuation part independent of the boundary conditions and a mean field part depending on 
the boundary conditions (which may fluctuate themselves, see the cases (A), (B2), and (B3)): 

2 

-Zfiuc = Z 1 / Vv J S(v(x.i))cxp 



(Vu) 



u 

a! 



(2) 



2 

2 / j \~ int 
m>l 
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Z mi = n / T>U exp E _£„ ^A^) ( 9 ^mf(xi))| exp|- ^0^ j .(5) 

The first exponential in Z m f stems from applying Gauss' theorem to the capillary wave Hamil- 
tonian Ti. C w[u m f + v]. In this term d n u m f denotes the normal derivative of the mean-field 
solution towards the interior of the circle dSf, and dli is the infinitesimal line segment on 
dSf. In the following we provide only the main steps of the subsequent evaluation of both 
the fluctuation and the mean-field part. For more details the reader is referred to Ref. [14]. 

Fluctuation part. The fluctuation part appears in all cases introduced above. In 
the case (Bl) it constitutes the full result for the partition function because in that case 
« m f = and Z m { = 1 . The 5- functions in the fluctuation part of the partition function can 
be removed by using their integral representation via auxiliary fields ipii^-i) defined on the 
interface boundaries dS® [13]. This enables us to integrate out the field u leading to 



-Z fluc = /ft Vtpi exp < / Mi (f d£ j i(}i(xi)G(\x i -x j \)t() j (x j )\ . (6) 

J ~r z i /rf, Jos Jos 



i=l 

Here, we introduced Green's function G(x) = A' (|x|/A c )/(27r) of the operator (—A + X~ 2 ) 
where Kq is the modified Bessel function of the second kind. In this form, the fluctuation 
part resembles 2d screened electrostatics: it is the partition function of a system of fluctuating 
charge densities %pi residing on the contact circles. For large d/ro it can be calculated by a 
multipole expansion [10]. In this limit and for X c /d ^oowc find for the fluctuation force: 

F Ruc = k B T^-\nZ auc - / +0(d- 3 ), -»1, ^--0. (7) 

od 2 d\n(d/r ) r A c 

Note, however, that here the limit A c — * oo is attained slowly with a leading correction term 
of the order l/(dln(A c /ro)), and that in this limit the free energy difference T{d) — J-(oo) 
is actually ill-defined (J-(d) ~ lnln(d/ro)) and therefore the effective colloidal interaction in 
case (Bl) - fixed colloids and pinned interface - is only meaningful for a finite capillary length 
A c [15], similar to a free two-dimensional interface the width of which is determined by the 
capillary wave fluctuations and diverges logarithmically ~ A c . 

In the opposite limit of small surface-to-surface distance h = d — 2tq the fluctuation 
force can be calculated by using the well-known result for the fluctuation force per length 
hd{h) = — ksT tt / (12h 2 ) between two lines a mean distance h apart [13], together with the 
Derjaguin (or proximity) approximation [16]: 

7rfc B T f ra 1 ro/h^oo ir 2 r)J 2 _ 1/2 



(h + 2r -2^r 2 -y 2 y 



This strong increase as h — > is a consequence of the finite (mesoscopic) size of the colloids 
and is missed if the colloids are approximated as pointlike objects [6]. It also turns out (see 
below) that this increase is not cancelled by any mean-field contributions. Thus it constitutes 
the dominant short-ranged contribution in all cases (A) and (B1)-(B3). For intermediate 
distances d between these limits -Ffluc must be evaluated numerically. In eq. © the integral 
over the auxiliary fields ipi can be carried out because they appear only quadratically in 
the exponent. The resulting determinant is divergent and requires regularisation. However, 
the derivative of its logarithm with respect to d (corresponding to the force) is finite and 
convergent in a numerical analysis (see Refs. [3, 14] for further details). Actual numerical 
results, which recover the limiting behaviors given by eqs. and JHJ, are presented in FigE] 
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Mean-field part. - The calculation of Z m f (Eq. JSJ) requires to determine the solution 
of the differential equation (—A + A~ 2 )u m f = for the (fluctuating) boundary conditions 
Mmf(xi) = /i(xj) for x, £ dS® and u m f(x) — > for |x| — > oo. We use the superposition 
ansatz u m { — ui + ui where m = J2 m K m {ri/\ c ) [Ai m cosmipi + Bi m sinm^i] is the general 
mean-field solution in R 2 \ S°. The solution has to match to the boundary conditions at 
both circles dSf and dS®- This can be achieved by a projection of U2 onto the complete set 
of functions on dS®, {cosmy>i,sinmy>i}, and vice versa. Equating this expansion with the 
multipole moments fi = (Pio, Pn,Qn, ■ ■ ■) of the values /i(xj) of the field u(x) at the contact 
lines leads to a system of linear equations for the expansion coefficients {A im , B im }. This 
system can be solved analytically within a systematic 1 jd expansion or numerically, observing 
rapid convergence (even for small d), allowing us to truncate the expansions at m max « 20. 
The mean field part of the partition function can then be written as 



cx P j-^_> 



^(Pio-hA 



where 7i[u m f] is a symmetric quadratic form of the multipole moments of the contact lines, 



T 



so that the d-dependent part of Z m { is given by det E. Before discussing the analytic structure 
of the dependence on d, we recall that the integration measure Vfi differs for the cases (A), 
(B2) , and (B3) . In all cases, the large d expansion of the mean-field part of the partition 
function (with A c — ► oo) leads to a repulsive effective force between the colloids: 

Due to Z = Zfi uc Z m f the total effective force is F — Fa uc + F m {. The leading terms in _F m f and 
-Ffluc (Eqs. (|llfl and Q) cancel. The same holds for the first subleading terms, and in the cases 
(A) and (B3) also for the second-next subleading terms. Due to these cancellations it turns 
out that in all four cases for large distances the effective force is attractive with the leading 
term stems from the fluctuation part. This can be summarized as follows (— 3> 1 , ■£ ► 0): 



F -> - 



8akBT fro\ 9 J (A) no pinning - fluctuating contact line: a = 1 



/'o \ d J \ (B3) pinning - colloidal height and tilt fluctuations: a = 9 



F — > — (-7) , (B2) pinning - colloidal height fluctuations 



ro V d 

p _ » — -j— — - , (Bl) pinning -fixed colloids. (12) 
2 dm{d/ro) 

In the opposite limit, h = d — 2R — > the effective mean-field force F m f is repulsive. In case 
(B2) the leading behaviour for h — > can be estimated analytically [14]: 

FnrfCft-O)**^ (13) 

We find numerically that the inclusion of higher (beyond the zeroth order) multipole moments 
of the contact line, as relevant for the cases (A) and (B3), slightly changes the prefactor in 
Eq. I|13|l but does not alter the l/h behaviour. The increase of the mean-field force for h — > 
is weaker than that of Ffi uc so that for all cases the leading behaviour of the total effective 
force at small distances is also attractive and given by Eq. (JSJ. 
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Fig. 1 - (a) Numerically calculated thermal total Casimir force F (symbols) for various boundary 
conditions compared with the analytic results in eq. 1121 (lines). The inset shows the results for the 
two contributions to the short range regime: numerically calculated fluctuation contribution (— Fn U c) 
and mean-field contribution (-F m f) compared with the analytic expressions in eqs. © (dashed) and 
11311 (dotted), respectively, (b) Sketch of the various boundary conditions. For (B1)-(B3), the disks 
may be replaced by Janus spheres. 

Discussion. The restrictions which two rotational symmetric colloids trapped at a 
fluid interface impose on its thermally excited capillary waves, lead to a thermal Casimir force 
between them. At large separations d 3> r$ there is an interesting interplay between the 
attractive interaction due to the interface fluctuations and a repulsive interaction generated 
by the mean- field contribution to the fluctuating boundary conditions. This results in a 
cancellation of leading terms up to a certain order in which depends on the specific 
type of boundary conditions. In the opposite limit of small separations h = d — 2tq <C t"o, 
the effect of the boundary conditions is much less pronounced, and the resulting force is 
dominated by -Fh uc ~ ft, -3 / 2 (attractive) compared with F m f ~ h" 1 (repulsive), leading to 
a strong attractive Casimir interaction in this regime (see Fig. Note that for colloid- 
colloid separations h of the order of the molecular length scale a of the fluids, the capillary- 
wave model is no longer valid. In this limit the total effective force stays finite and can be 
understood by taking into account the formation of a wetting film around the colloids [17]. 
At small separations typically also van der Waals (vdW) forces become important, leading 
to a strong tendency of colloid aggregation if not compensated by a repulsive interaction. 
For spherical colloids at small separations and 6 = tt/2 the well-known Derjaguin result for 
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the vdW force is given by F^^ = — Ah/(127t) R/h 2 where the (effective) Hamaker constant 
An is determined by the frequency dependent dielectric permittivities of both the colloids 
and the two fluid phases. shows an even stronger increase for h — ► and hence will 

dominate the Casimir force (JBJ since the Hamaker constant is typically An rj 1 . . . lOfceT. 
However, refractive index matching between colloids and fluids can result in a significantly 
smaller value of A H . For An/k B T < (7r 3 /4) y/h/R we have |F fluc | > [Ff^l and therefore 
in this distance regime an increased influence of the fluctuation induced force on particle 
aggregation. In the case of two cylindrical disks with height H the Derjaguin approximation 

leads to F^w {H/h 2 ) (R/h) 1 / 2 if h < R, H. However, for very fliin disks (H < /i « R, 

c.f. cases (B1)-(B3)), integration over atomic pair potentials ~ ?'~ 6 , which give rise to the 
vdW force, results in 



For these systems F vdw /F fluc = 15(An/k B T)(H/h) 2 < 1, so that for H < h < R the 
driving mechanism for flocculation is entirely given by the fluctuation induced force. Thus, 
the fluctuation induced force can strongly enhance the tendency of colloids at interfaces to 
flocculate. This effect is independent of material parameters, as long as the system is in the 
capillary wave regime, i.e., for h 3> o. In order to discriminate the fluctuation induced force 
from the vdW force we propose three experimental scenarios. The first one is to weaken 
the vdW force by lowering the Hamaker constant via index matching which increases the 
importance of the fluctuation force. Second, for thin disks, the fluctuation induced force 
dominates, independent of the dielectric properties of the materials. A third mechanism 
consists of reducing the vdW force by a light polymer coating of the colloids so that in eq. (JSJ 
h = h sts is given by the surface-to-surface distance between the polymer shells, whereas the 
distance h = h COIC > h sts between the cores of the two colloids enters into . Thus for a 

sufficiently thick polymer shell the fluctuation induced force will be the dominant interaction 
at small distances. 
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